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INTRODUCTION 


Continuous  wave  (cw)  chemical  lasers  generally  operate  at  pressures  of 
the  order  of  1  Torr  in  order  to  achieve  good  efficiency.  At  this  pressure 
level,  the  gain  medium  is  inhomogeneously  broadened.  In  addition,  the  lasing 
process  tends  to  modify  the  gain  medium  so  that  the  lasing  particles  are 
neither  in  translational  nor  rotational  equilibrium.  Nonetheless,  early 
analytic  models  of  cw  chemical  laser  performance  assumed  that  the  lasing 
medium  was  in  translational  and  rotational  equilibrium  (e.g.,  Refs.  1  and 
2).  In  these  models,  reasonable  estimates  were  provided  for  net  output  power, 
but  the  6pectral  content  was  not  predicted.  The  latter  requires  consideration 
of  finite  translational  and  rotational  relaxation  rates. 

Subsequently,  analytic  models  and  numerical  codes  were  developed  that 

3-5 

included  either  rotational  nonequilibrium  or  translational  nonequilib- 
— ft 

rium.  A  recent  review  of  rotational  nonequilibrium  rate  processes  and 

models  is  given  in  Ref.  9. 

The  combined  effect  of  translational  and  rotational  nonequilibrium  has 

received  less  attention.  An  analytic  model  for  a  low  pressure  CC>2  laser  is 

presented  in  Ref.  10,  whereas  a  mode2  for  a  generic  molecular  laser  with 

applications  to  CO2  and  cw  chemical  lasers  is  presented  in  Ref.  11.  A 

numerical  code  that  treats  both  translational  and  rotational  nonequilibrium 

12 

has  been  developed  by  D.  Bullock  and  co-workers;  limited  numerical  results 

13 

have  been  published. 

6” ft 

The  present  study  generalizes  a  previous  simple  model  in  order  to 
include  rotational  as  well  as  translational  nonequilibrium  effects.  The 
object  is  to  provide  analytic  expressions  for  cw  chemical  laser  amplifier  and 
oscillator  performance  as  well  as  to  delineate  the  parameters  that  charac¬ 
terize  nonequilibrium  effects.  The  appropriate  equations  are  first  deduced. 
Equilibrium  and  nonequilibrium  solutions  are  then  obtained.  Symbols  are 
defined  in  Appendix  A. 
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II.  FORMULATION 

Equations  are  deduced  that  define  effects  of  translational  and  rotational 
nonequilibrium  on  cw  chemical  laser  performance. 

A.  Flow  model 

A  cw  chemical  laser  is  illustrated  in  Fig.  la.  The  present  simplified 
mixing  model  is  illustrated  in  Fig.  lb.  The  reactants  are  assumed  to  be 
premixed  but  do  not  react  until  a  flame  sheet,  y^Cx),  is  reached.  The  flame 
sheet  shape  is  specified,  a  priori,  from  diffusion  theory.  The  streamwise 
station  where  the  flame  sheet  reaches  the  channel  center  line  is  denoted  x^, 
and  characterizes  the  diffusion  rate.  The  width  per  semichannel  and  the 
number  of  semichannels  are  denoted  w  and  ng„,  respectively.  Laser  radiation 
is  in  the  +y  direction. 

B.  Distribution  functions 

Let  nv(J,v)dv  denote  the  number  of  particles  (moles/craJ)  with  vibrational 
energy  level  v  and  rotational  energy  quantum  number  J  that  are  resonant  with 
radiation  of  frequency  v.  The  following  notation  is  introduced 


/  nv(J,v)dv  =  ny(J) 


v  n  ( J )  =  n 
J  v  v 


[  11  =  nT 

L  v  T 


For  radiation  in  the  ±y  direction,  the  resonant  frequency  v  is  related  to 
particle  thermal  velocity  v^  by  the  Doppler  relation 

1  -  (2) 
vc 
o 

where  v  is  the  resonant  frequency  for  particles  at  rest.  We  neglect  the 

dependence  of  v  on  v,J. 

o 

For  particles  with  a  Maxwellian  thermal  velocity  distribution 


-  .  nv(J*u)  [  ( 4  L  n  2  )  /TT  ) 


n  ( J) 

V 


exp  -(4tn2)^ 


X3  w  jutx  *  *  *  k  w~a  'wcd.xu'r*  iruxv  *v  vv»vwv^.rWviY.AAv,JV‘  v<  ’/  ^^jrmjrmjr*3r*jr^jr«ji 


where  Av  is  the  characteristic  Doppler  width  (Appendix  B).  The  value  of  p  at 
d 

line  center  is  denoted  pQ  and  is  used  to  characterize  the  reciprocal  of  the 
Doppler  width.  Thus 


p  =  [  (4Zn2)/x  ]  ly,^/Av 
o  d 


We  define  p  =  p/p  and  note 
r\>  r  o 


00  , 
/  P  dv  =  ~ 


We  now  consider  the  case  of  rotational  equilibrium  and  neglect  the  effect  of 
vibrational  energy  level  on  the  characteristic  rotational  energy  temperature 
Tr.  The  fraction  of  particles  in  rotational  energy  level  J  is  denoted  fj  and 
is  found  from  (Appendix  B) 

_  nv(J)  ( 2 J  +  l)exp[-J(J  +  1)9r] 

fJ  =  =  Z(2J  +  1  )exp [ -J ( J  +  l)a  1  (6) 

v  J  R 

where  the  summation  is  from  J  =  0  to  J  *  00 .  Let  f  denote  a  reference  value 

r 

of  f,  and  introduce  f  ,  =  f , 'f  .  It  follows  that 
J  J  J  r 

I  f  -  1  ;  r  f  -  1/f  (7) 

J  J  J  J  r 

Convenient  values  for  f  are  noted  in  Appendix  B.  Finallv,  we  observe  that 
the  gain  at  frequency  v  for  laser  transition  v,J  can  be  expressed 

gv,J(v)E%,J  '  (nv(J.v'))u  -  [nv(J,v')]£|  L(v  -  v')dv>  (8a) 

where  subscript  v,J  refers  to  lower  laser  level  values,  subscripts  u 

and  l  denote  upper  and  lower  laser  level  values,  respectively,  and 

L(v  -  V')  =  rLl  +  1  (8b) 

h 

where  Av  is  the  homogeneous  line  width.  The  evaluation  of  Av,  and  the  cross 
h  h 

section  0  is  discussed  in  Appendix  B.  Numerical  estimates  are  provided  in 
v,  J 

Table  I. 


^  -C  -  JA.  7 LaVa  VWVaV  Vm  mV  mV  mV  A 


J  m"  U*  «<*'-*•* < 


C.  Laser  equations 

In  order  to  simplify  the  mathematical  development,  we  assume  a  ”Q"  type 
laser  transition 


v  +  1 ,  J  v,  J 


rather  than  the  “P"  type  laser  transition 


v  +  1,  J  *  1  +  v,J 


appropriate  for  cw  chemical  lasers.  This  approximation  is  consistent  with 
other  simplifying  assumptions  used  in  the  present  model.  In  addition,  we 
consider  a  two-vibrational  energy  level  model  and  denote  the  upper  and  lower 
levels  by  subscripts  2  and  1,  respectively. 

The  variation  of  upper  level  number  density  with  streamwise  distance,  in 
the  present  model,  is  found  from 


u  d[n2(J,v)yf; 


'  'j  J,Vf>  ,  ,,  , 

*  u  TT  d>""~  kcd"2<J,V) 


+  ktr[pn2(J)  -  n2 ( J ,v ) ]  +  krrp[fjn2  -  n2(J)] 


(10a) 


-  — 3 —  [n?(J,v)  -  n,  ( J ,v  )  ]  l  L(v  -  v.  ,)I.  , 
v,J  Z  1  j  J’3  J,j 

where  e  is  the  energy  per  mole  of  photons.  The  terms  on  the  right-hand 

v ,  j 

side  of  the  equation  represent  the  effects  of  chemical  pumping,  collisional 
deactivation,  translational  cross  relaxation  (assumed  to  be  proportional  to 
departure  from  translational  equilibrium),  rotational  cross  relaxation 
(assumed  to  be  proportional  to  departure  from  rotational  equilibrium),  and 
stimulated  emission  and  absorption.  For  the  term  involving  ktr,  it  is  assumed 
that  within  a  given  rotational  level  J,  the  creation  of  n2(J,v)  is  propor¬ 
tional  to  the  departure  from  translational  equilibrium.  For  the  term  involv¬ 
ing  krr,  it  is  assumed  that  rotational  relaxation  into  level  J  results  in 
particles  with  a  Maxwellian  velocity  distribution.  Equation  (10a)  provides 
the  proper  limits  n2(J,v)  +  pn2(J)  and  n2(J)  +  ^jn2  as  ^  *  “  ar|d  +  "> 


respectively.  It  is  assumed  that  the  chemical  reaction  creates  only  upper 

level  particles  and  that  these  are  in  translational  and  rotational 

equilibrium.  The  notation  I.  .  and  v.  refers  to  the  resonator  mode  for 

J » J  d, j 

the  v  +  1,  J  +  v,J  transition.  The  corresponding  expression  for  the  lower 
laser  level  is 


u  dln^J.vtyj] 


kcd  p  fjn2  +  k  [pn^ (J)  -  n^J.v)] 


+  krrp[fJn1  -  n^J)] 


(10b) 


+  [n2(J,v)  -  n^(J,v  ) ]  l  L(v  -  Vj 

v,J  j  * 2  ,;1 


The  first  term  on  the  right-hand  side  follows  from  the  assumption  that  colli- 
sional  deactivation  of  n2  results  in  an  equilibrium  distribution  of  n^(J,v). 
The  following  nondimensional  quantities  are  introduced 


n2(J,v)yf 


n  f  p  w 
r  r  o 

n;  =  N__  ±  N, , 

Jv  2Jv  lJv 


N“  =  N2  ±  N 


R  -  k  /k  . 
t  tr  cd 


n  f  p  ivLo  ,w 
r  r  o  h  v,J 


n2(J)yf 

(11a) 

n  f  w 
r  r 

n2yf 

n  w 
r 

(lib) 

N2J  *  N1J 

(11c) 

IJ,j/IS,L 

(lid) 

k  /k 
rr  cd 

(lie) 

W 

(Ilf) 

/*  N L(v  -  „•>  p- 
J  Jv  Av, 

-oo  h 


S  •  tt  p  f  Av, 

ro  r  h  t 


Olg) 


where  nr  is  a  reference  value  of  n2  and 


13 


i  . ' |  g'f  #*•  #’4  W*a  1^1  a 


.**  - .ta.U  JB  .fa*. 


•a* Ji*  a*  A? :  v  i 


Ic  -  e  Tk  /(2o  .)  (12) 

S,L  v,J  tr  v,J 

is  an  intensity  that  characterizes  line  shape  distortion  resulting  from  satu¬ 
ration. 

Substitution  of  normalized  variables  into  the  sum  of  Gqs.  (10a)  and  (10b) 
and  subsequent  integration  with  respect  to  v  and  summation  with  respect  to  J 


indicates 


dNJv  dNT  +  + 

~~  -  p  f  T  — —  +p  f,N,  -N„  +  R  (p  N  -  N,  ) 

d?  Kv  J  dc  Kv  J  2  2Jv  t  y\>  J  Jv 


+  Rr<\,<fjN+  -  NJ> 


dNJ  dNT  +  + 

“dT  ‘  f J  IT  *  «jn2  -  N2J>  +  W  -  NJ  > 


N*  -  ». 


The  difference  between  Eqs.  (10a)  and  (10b)  indicates 

dNJv  dNT  -  - 

-  p  f.  --  ~  (P  f  TN«  +N„  )  +  R  (p  N.  -  N.  ) 
d£  v  J  d;  'v  J  2  2Jv  t  rv  J  Jv 


+  Vv(fJN"-V  -RtNJv  l  L(V 

dNj  dNT 

-dT=f J—  (fJN2+N2J)  +Rr(fJN  -MJ> 


(13a) 


(13b) 


(13c) 


( 14a) 


p  4v  R  y  G.  .1.  . 

H°  h  C  j  J'l 


dN_  J  - 

-  (N_  +  N  )  -  -  7  £  G .  ,1  .  , 
d;  d;  T  v  Lj  ^  J,J  J,j 


(14b) 


(14c) 


Equations  13  and  14  can  be  solved  for  Njv i  Nj>  and  N”,  if  the  chemical  pumping 
rate  dNT/d?  is  specified  and  if  Gj  j  (oscillator)  or  I  j  j  (amplifier)  is 
specified.  For  laminar  mixing 


nt  -  (;/;D) 


mm- 


I 


1 


where  -  k  ,x_/u  is  the  normalized  diffusion  distance.  For  a  Fabry-Perot 

D  Cd  D 

resonator  where  each  mirror  has  a  reflectivity  Rm,  gj  jyfngc  *  -inR^,  and 


■  constant 


n  n  f  p  Av.c  w 
sc  r  r  o  h  r 


The  output  power  per  semichannel,  released  up  to  station  x,  is  denoted  P  and 


is  found  from 


- - —  ■  oT  l  l  T/e  )  /  G.  .1.  .dc 

nrwuer  2^  j  i  V»J  r  0  J,:I  J,j 


(17a) 


If  we  assume  e  ,  =  e  ,  Eq.  (14c)  indicates 
v,J  r 


2P  -  N  -  N  -  /  (N  +  N“)d? 


The  solution  of  these  equations  is  discussed  in  the  following  sections. 


(17b) 


'  v„v'N\ 


■I 


H*  I.- 


III.  EQUILIBRIUM  CASES 

The  assumption  of  translational  and  rotational  equilibrium  corresponds  to 
taking  the  limits  Rt  +  ®  and  Rr  -*•  00 ,  respectively,  in  Eqs.  (13)  and  ( 1 A ) .  The 
results  are 


N't  N't 

-4-  -  /  -  N+ 

pvfJ  fJ 


(18a) 


I 

I 

I 

£ 


N.  N. 

Jv  J 

”  7—  “  N 

p  f,  f  . 

*v  J  J 


as  expected  and 


sitN'>  -  fn  j 


Also,  from  Eqs.  (18b)  and  (Ilf),  in  the  limit  Av^  «  Av^, 


2  GJ,-1 


exp[  -4in2( 


v .  -  v  -  Av. 

J  Oi  A  r  ,  .  «  f  h 


W7V)  [fo(jJ)] 

d  a 


(18b) 


(19a) 


(19b) 


In  the  case  of  an  oscillator,  threshold  gain  is  specified  (Gj^j  *  Gc),  and 
Eqs.  (19a)  and  (19b)  are  solved  for  Ij  j.  In  the  present  two-vibrational 
level  model,  one  lasing  transition  J,j  occurs,  namely  the  transition  with 
highest  gain.  (In  a  multivibrational  level  model,  there  is  one  lasing  transi¬ 
tion  for  each  upper  vibrational  level.)  Amplifier  solutions  are  obtained  by 
specifying  I.  .  in  Eq.  (19a).  Equations  (19a)  and  (19b)  correspond  to  the 
equilibrium  model  presented  in  Ref.  1. 


m 


$ 

1 

i>J»| 

w 

Ml 

ts 

I 

K? 


‘•*2 

$ 

A 

A 
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IV.  NONEQUILIBRIUM  CASES 


We  now  consider  effects  of  translational  and  rotational  nonequilibrium. 
The  quantities  Rr  and  Rt  are  large  in  cw  chemical  lasers  (Table  I  and  Ref. 
9).  Therefore,  we  consider  the  limit 

Rr  »  1,  Rt  »  1  (20) 

A  similar  approximation  was  Introduced  in  Ref.  11.  Simplified  laser  equations 
are  deduced.  Amplifier  and  oscillator  solutions  are  then  obtained. 

A.  Simplified  equations 


1.  Limit  Rr  »  1,  Rt  »  1 

Recall  that  number  density  variables  in  Eqs.  (13)  and  (14)  have  been 
normalized  to  be  of  order  1.  If  terms  of  order  Rr  *  and  Rt~*  are  neglected, 
Eqs.  (13a)-(13c)  indicate 


Thus,  N*  and  N*  retain  translational  and  rotational  equilibrium  distributions 
Jv  J 

in  the  present  approximation.  Similarly,  Eqs.  (14a)-(14c)  become 
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Convection,  chemical  pumping,  and  collisional  deactivation  terms  no  longer 
appear  in  Eqs.  (22a)  and  (22b).  In  the  present  limit,  the  rate  of  increase  of 
(or  Nj)  resulting  from  cross  relaxation  is  just  equal  to  the  rate  of  loss 
of  N~  (or  Nt)  resulting  from  stimulated  emission.  Thus  the  cross  relaxation 
and  stimulated  emission  processes  are  in  equilibrium,  and  other  rate  processes 


are  too  slow  to  affect  this  equilibrium.  The  convection,  chemical  pumping, 
and  collisional  deactivation  terms  are  retained  in  Eq.  (22c),  which  represents 
global  conservation.  Substitution  of  Eq.  (22a)  into  Eq.  (Ilf)  indicates 
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Equation  (23)  exhibits  the  variation  of  gain  with  frequency  for  a  given  tran¬ 
sition  v,J.  The  quantity  Kj(v')  is  the  normalized  line  shape  and  is  seen  to 
depend  only  on  Av  Au , ,  and  IT  The  coefficient  of  K  (v  ' )  in  Eq.  (23a)  is 
an  amplitude  function  and  requires  a  knowledge  of  the  number  densities  N-  and 
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Evaluation  of  the  line  shape  can  be  simplified  by  introduction  of  the 
variables 
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For  cases  where  Av,  <<  Av  ,  the  value  of  K.(v)  at  v.  is  affected  only  by  the 
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laser  intensity  at  v  and  is  found  from 
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where  D(Xj)  is  the  Dawson  integral.^  The  quantity  X^  can  be  replaced  by  X  in 
Eqs.  (25a)  and  (25b)  when  I.  .  ■  0.  In  this  case,  Eqs.  (25a)-(25c)  provide 

J  >  J 

expressions  for  the  line  shape  Kj(v). 

2.  Limit  Rr»l,  Rt»l,  |  Rr /Rt  “  l|  (Avh/Avd)  «  1 

Equations  (22a)-(22c)  can  be  further  simplified  if  it  is  assumed  that 
Rr>>! ,  Rt»l  and 
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Equations  (22a),  (22c),  and  (23a)  become 
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These  equations  can  be  evaluated  without  consideration  of  Nj  ,  which  is  found 
f  rom 
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Note  that  Nj~  is  reduced  below  its  equilibrium  value  by  an  amount  which, 
for  Kj(Vj)  -  1 ,  is  proportional  to  IjAv^/Av^.  Equations  (27a)-(27c)  are  the 
same  equations  that  result  when  rotational  equilibrium  is  assumed. 

Equations  (20)  and  (26)  are  realistic  for  low  pressure  cw  chemical 
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lasers,  because  R  ~  R  »  1  and  Av,/Av  ~10  p(Torr)  <<  1  in  these  lasers 
t  r  h  d 

(Table  I).  With  an  increase  in  pressure,  Eq.  (26)  remains  valid  if  it  is 


1  Is  consistent  with  estimates 


assumed  that  R  /R^  m  1.  The  assumption  R  /R^  * 
r  t  r  t 

in  the  range  0.5  <.  Rr/Rt  <_  10  in  Refs.  9  and  11  and  with  the  simplified  nature 
of  the  present  model.  The  present  results  suggest  that  a  first  estimate  for 
the  performance  of  cw  chemical  lasers  can  be  obtained  by  assuming 
translational  nonequilibrium  and  rotational  equilibrium  [Eqs.  (27a)-(27c)] . 
The  quantity  Nj~  is  then  found  from  Eq.  (27e). 


B.  Amplifier 


We  consider  a  multiline  amplifier  (Fig.  lb)  in  the  limit  given  by  Eqs. 
(20)  and  (26)  and  assume  that  the  input  intensity  of  each  transition  Ij  j  is 
specified  and  is  independent  of  £ .  The  quantity  B  in  Eq .  (27)  is  then  a 
constant.  For  laminar  mixing,  integration  of  Eq .  (27b),  together  with  Eq. 
(27c),  yield 
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where  D(  )  is  again  the  Dawson  integral.  Equation  (28)  provides  the  varia¬ 
tion  of  N~  with  streamwise  distance  for  various  values  of  the  parameter  B. 
Numerical  results  are  plotted  in  Fig.  2.  The  maximum  points  on  these  curves 
are  denoted  by  subscript  m  and  are  included  in  Table  II. 

Equation  (28)  also  provides  the  variation  of  gain  with  frequency,  stream- 
wise  distance,  and  saturation.  The  line  shape  corresponding  to  a  single  laser 
transition  at  frequency  Xj  =  0  and  Xj  »  0.40  is  shown  in  Figs.  3a  and  3b, 
respectively.  The  ratio  of  homogeneous  to  inhomogeneous  broadening  is  assumed 

to  be  Av/Av  =  0.024,  and  the  corresponding  low  saturation  hole  size, 

1  /  2  11  “ 

(£n2)  Av./Av,,  is  indicated.  Hole  size  depends  on  Y.  ,  [Eqs.  (25a )-(25c ) ] , 
and  the  significant  increase  of  hole  size  with  saturation  is  evident.  The 
variation  of  local  gain  with  frequency  at  5  r  0.305,  for  a  single  laser 
transition  of  Xj  =  0.40,  is  shown  in  Figs.  4a  and  4b  for  Av^/Av^  *  0.024  and 
0.048,  respectively.  The  station  c  *  0.305  corresponds  to  the  6treamwise 
location  where  the  zero  power  gain  is  a  maximum  (Table  II).  Increased 
saturation  is  seen  to  depress  the  entire  gain  curve.  This  is  a  consequence  of 
cross  relaxation. 
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STREAMWISE  DISTANCE,  f  =  kc.  x  u 

Fig.  2.  Variation  of  Inversion  Number  Density  with  Streamwise  Distance  for  Case 
of  Uniformly  Illuminated  Amplifier  with  Laminar  Mixing.  Parameter  B  is  a 
measure  of  saturation  level. 
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TABLE  II.  Maxima  for  amplifier  with  uniform  incident 
radiation  and  laminar  diffusion.3 
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Fig.  4.  Variation  of  Gain  with  Frequency  at  Streamwise  Station  £  «  0.305  in  a 
Uniformly  Illuminated  Amplifier  with  a  Single  Transition  at  X,  *  0.4. 
Laminar  mixing  and  homogeneous  widths  Av,  ,/Av  -  0.024,  0.048  are  assumed 
[Eq.  (23a)];  (b)  case  Av  /Av  *  0.048. 
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Reference  15  reports  experimental  measurements  of  the  decrease  in  line 
center  gain  at  streamwise  station  x  ■  0.4  resulting  from  amplifier  radiation 
at  frequency  Xj  “  0.40.  Figure  3a  and  Table  Id  can  be  used  to  provide 
theoretical  estimates  for  line  center  gain  variation  with  amplifier  input 
intensity.  The  resultant  estimates  are  included  in  Fig.  5  and  indicate  good 
agreement  with  experiment.  [The  authors  of  Ref.  15  interpreted  the  present 
line  shape  estimates  (Fig.  3)  as  line  center  gain  estimates  (Fig.  4)  and 
incorrectly  concluded  that  the  present  model  does  not  properly  evaluate 
saturation  effects.] 

The  downstream  end  of  the  positive  gain  region  is  denoted  and  is  found 
by  equating  N~  to  zero.  The  result  is 
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The  net  output  power  is 


ld/2?  =  [B/(l  +  B)k  1/2[1  -  (2/3  )C  ]  (30) 
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Corresponding  values  of  B,  and  Pg  are  included  in  Table  II.  For  a  sat¬ 

urated  (B  00 )  amplifier 
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Equations  (28)— (31)  are  identical  in  form  to  the  corresponding  results  pre¬ 
sented  in  Ref.  1.  The  parameter  B  in  these  equations  replaces  the  parameter 
I<2  in  Ref.  1.  The  previous  results  neglect  hole  burning  effects  on  Kj  j. 

The  single-line  amplifier  solution  provides  a  convenient  basis  for 

investigating  the  saturation  process  in  cw  chemical  lasers.  The  effect  of 

saturation  on  normalized  line  shape  is  determined  by  the  parameter  for 

an  inhomogeneous  medium,  the  intensity  Ij  j  *  1  results  in  a  reduction  of  Kj  j 

by  a  factor  of  2"*^  [Eqs.  (24d)  and  (25b)].  The  effect  of  saturation  on 

power  extraction  is  characterized  by  the  parameter  B,  as  indicated  in  Table 

II.  In  the  limit  B  -*■  ®,  all  available  power  is  extracted  by  a  single  line. 

For  the  case  of  a  single  lasing  transition,  a  power  extraction  saturation 

intensity  I  can  be  defined  by  Eq.  (27d) 
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5.  Effect  of  Input  Intensity  on  Line  Center  (v  =  vQ)  Gain  at  Fixed 
Streamwise  Station  for  Uniformly  Illuminated  Amplifier.  Circles  de¬ 
note  experimental  data  from  Ref.  15  that  were  taken  at  x  =  0 . ^  cm 
with  a  gain  medium  at  p  =  6  Torn  and  a  laser  transition  corresponding 
to  v,J  =  1,8  and  Xj  =  0.4.  The  symbol  f  denotes  present  estimate  for 
gain  based  on  line  center  values  from  Fig.  4a  and  L  =  46  W/cm'T 
The  estimate  for  L  ^  is  obtained  by  interpolation  from  Table  la  by 
recalling  l  -  p2  and  by  assuming  p  =  6  Torr  and  T  =  900  K.  Note 

also,  Table  id  indicates  =  0.025,  which  is  consistent  with 

use  of  Fig.  4a. 
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where 


0  “  2B/(Sf/jjKj>j) 


Values  of  Ic  p,  corresponding  to  B  ■  1.5  and  B  ■  7.0,  are  included  in  Table  I 

and  are  seen  to  be  large  when  S  is  small.  For  cases  with  N  strong  laser 

transitions  and  Av.  «  Av  ,  the  power  extraction  saturation  intensity  is 

n  c 

approximated  by  Eq.  (32)  with  0  replaced  by  0/N. 


C.  Oscillator 

We  now  consider  a  Fabry-Perot  oscillator  in  the  limit  given  by  Eqs.  (20) 

and  (26)  and  further  assume  that  Av,  «  Av  and  Y,  .  <<.  1.  These  assumptions 

n  C  J  *  J 

simplify  the  line  shape  [Eq.  25b)].  We  also  assume  that,  within  each  v,J 
lasing  band,  one  laser  frequency  corresponds  to  v^;  the  other  laser  frequen¬ 
cies  are  then  symmetric  about  v  ,  and  the  present  summation  with  respect  to  j ,, 

®  t 

rather  than  thermal  velocity,  needs  no  modification.  The  intensity  Ij  j  then 

represents  the  sum  of  the  intensities  in  the  +y  and  -y  directions.  The 

case  Av  <<  Av,  is  treated  in  Ref.  7. 
c  h 

The  threshold  gain  is  denoted  Gc  and  is  a  constant.  The  gain  for  each 
lasing  transition  is  then,  assuming  o  ,  -  o 
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Equations  (25b)  and  (33)  yield 
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where  N”  is  obtained  from  Eq.  (22c),  namely 
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Equation  (34a)  provides  the  dependence  of  Ij  j  on  fj  and  Xj  at  each  streamwise 
station.  The  number  of  rotational  levels  and  longitudinal  modes  that  reach 
threshold  at  each  streamwise  station  is  found  from  the  requirement  that  lj  j  > 
0  in  Eq.  (34a),  namely 
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The  number  of  lasing  transitions  increases  as  Gc/N~  decreases.  Equation  (34b) 
cannot,  in  general,  be  integrated  in  closed  form.  Two  subcases  are  treated:  a 
saturated  multiline  oscillator  and  a  partially  saturated  single  line  oscilla- 


1.  Saturated  multiline  oscillator 

Assume  that  Gc  <<  1  and  S  *  0(1).  In  this  limit,  N-  <<  1  and  N-/Gc^^  » 
0(1).  The  laser  is  saturated,  and  there  is  a  large  number  of  laser  transi¬ 
tions.  Equations  (34a)  and  (34b)  become 
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Corresponding  number  densities  are 
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For  cases  where  the  rotational  energy  levels  and  longitudinal  mode  frequencies 
are  closely  spaced,  the  summations  in  Eqs.  (36a)-(36d)  can  be  replaced  by 
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Equations  (36)  and  (26b)  provide  the  net  inversion  N“  and  the  lasing  intensity 
Ij  j  at  each  streamwise  station.  The  quantities 
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are  plotted  in  Figs.  6a  and  bb  for  Gc  -*■  0.  Longitudinal  mode  intensity  Ij  j 

is  inversely  proportional  to  Av  ,/Av  ,  which  is  a  measure  of  the  number  of 

d  c 

longitudinal  modes  with  a  significant  amount  of  power.  The  dependence  of  Ij 

on  9r  is  indicated  in  Fig.  6b.  The  number  of  active  lasing  modes  decreases  as 

Gc  increases.  Equation  (36c)  indicates  that  departs  from  the  equilibrium 

value  p  f  N  only  in  the  vicinity  of  each  lasing  frequency  v.  , .  The  depar- 
v  J  J >  j  __ 

ture  from  equilibrium  is  large  because  Ij  j  is  large.  The  departure  of  Nj~ 
from  the  equilibrium  value  f jN  is  of  order  l/tR^N  ). 

Because  the  laser  is  saturated,  the  output  power  and  mode  length 
1/2  1/2 

are  2'  P  -  2  /3  and  5  *  1/2,  respectively,  for  laminar  mixing.  Thus, 

D  e  e 

for  the  saturated  multiline  oscillator,  rotational  and  translational  nonequil¬ 
ibrium  impact  spectral  output  but  not  output  power. 

2.  Single  line  oscillator 

The  parameters  S  and  Gc  determine  the  degree  of  saturation.  In  order  to 
investigate  their  influence,  we  consider  a  case  where  line  selection  results 
in  a  single  laser  transition.  For  convenience,  consider  laminar  mixing,  fj  * 
1,  and  Xj  *  0.  Lasing  is  initiated  at  station  where  the  gain  first 

reaches  the  threshold  value.  Integration  of  Eq.  (34b)  with  IT  .  ■  0  indicates 
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ROTATIONAL  LEVEL  J 


Performance  of  Saturated  Multiline  Fabry-Perot  Oscillator  in  the  Limit  Rr 
=  R  »  1,  «  AX  c,  Yj  ,  «  1,  Gc  -  0,  and  S  =  0(1)  [Eqs.  (36)  and 

(37)].  (a)  Variation  Cof  longitudinal  mode  intensity  Ij  .  with  frequency 

X.  and  (b)  variation  of  intensity  I,  with  rotational  level  J. 
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which  agrees  with  Eq.  (28)  for  B  *  0.  Equation  (38)  provides  5  and  ?  N 
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for  a  given  value  of  G  .  Downstream  of  c  ^ 
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which  is  integrated  with  the  initial  condition  N~  -  N^-  at  c  -  Lasing 

terminates  at  the  downstream  station  where  N-  ■  N^~.  The  net  output  power 
is 

2;d1/2  p6  -  [c 1/2  -  (2/3)c 3/2]^e  -  /  6  ;d1/2n>  (*o) 

^  5  < 


Equations  (39)  and  (90)  have  been  evaluated  for  several  values  of  S  and 

1/2 

G  .  Output  power  is  given  in  Fig.  7  and  decreases  as  S  and  G  are 
D  C  c 

reduced. 

Oscillator  solutions  where  hole  burning  effects  are  neglected  correspond 
to  S  +  Ij  j  -*■  0,  and  SIj  j  *  finite.  The  number  density  in  the  lasing 
region  is  a  constant  given  by  N~  -  Ni~»  and  are  obtained  from  Eq. 

(38).  The  local  lasing  intensity  for  laminar  mixing  is 

9r  1/2  dP  _  1  1/2  l/2„- 

2CD  d^-~l/2  C  CD  Ni  (41) 

The  end  of  the  lasing  region  occurs  when  the  intensity  goes  to  zero,  or 
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The  net  output  power  is  found  from  Eq.  (90)  and  is  included  in  Fig.  7. 
Neglect  of  hole  burning  is  seen  to  overestimate  output  power. 


DISCUSSION 


The  parameters  that  characterize  laser  performance  are  liscussed  further, 
and  numerical  estimates  are  provided. 

The  parameters  introduced  in  the  present  study  may  be  expressed  in  the 


ktr  n2(J,v^ 

Kd  n2(J*v) 


Particle  Collision  Rate 
Particle  Deactivation  Rate 


(°v,J/cv,J)fJ,j  n2(J»v) 
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Resonant  Particle  Collision  Rate 
Net  Deactivation  Rate 
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J,j  Avnn2(J,v) 
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Net  Stlm.  Emission  Rate 
Net  Deactivation  Rate 


Net  Stlm.  Emission  Rate 
Net  Deactivation  Rate 
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(43c) 
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The  parameter  Rt  may  be  viewed  as  the  number  of  translational  collisions  an 

upper  level  particle  undergoes  before  it  is  col  1 isional ly  deactivated.  The 

parameter  Ij  j  represents  the  ratio  of  particle  stimulated  emission  rate  to 

particle  collisional  deactivation  rate.  Similarly,  S  represents  the  ratio  of 

the  collisional  deactivation  rate  of  particles  resonant  with  lj  j  to  the  net 

upper  level  particle  collisional  deactivation  rate.  Finally,  B  and 

G  S  £  l  I.  .  apply  to  amplifiers  and  oscillators,  respectively,  and  represent 
J  i 

the  ratio  of  net  stimulated  emission  rate  to  net  collisional  deactivation 
rate.  The  latter  ratio  characterizes  laser  output  power  and  efficiency. 

Numerical  estimates  for  parameters  and  rate  coefficients  are  included  in 
Table  I  for  a  fixed  stoichiometry,  300  <  T,  K  <  1200,  0.  132  <  p,  atm  <  13.2 
(0.10  <  p,  Torr  <  10.0),  and  v,J  ■  0,  6;  0,  9;  1,5,  and  1,  8.  The  pressure 
dependence  of  quantities  in  Table  I  is  given  by 


i,-*  \ 


and  permits  interpolation  of  data  in  Table  I  with  regard  to  pressure.  Values 
of  the  power  saturation  intensity  I  are  given  for  B  =  1.5  and  7.0.  These 

_  ^  9  * 

values  of  ^  correspond,  roughly,  to  achieving  one-half  saturated  power 
output  and  one-half  zero  power  inversion  number  density,  respectively.  Note 
that  Ig  p  is  relatively  insensitive  to  pressure,  in  contrast  to 


%  *  .V  A  \  .S  v 


VI .  CONCLUDING  REMARKS 


In  the  present  model  we  have  assumed  that  convection,  chemical  pumping, 
and  collisional  deactivation  rates  are  small  compared  with  translational 
relaxation,  rotational  relaxation,  and  stimulated  emission  rates.  Similar 
approximations  are  introduced  in  Ref.  11.  It  is  expected  that  corresponding 
simplifications  can  be  introduced  into  numerical  codes  in  order  to  facilitate 
solutions . 

The  further  assumption  |Rr/Rt  “  1 1 (Av h /Av ^ )  «  1  resulted  in  a  system  of 
equations  for  cw  chemical  laser  performance  that  are  independent  of  Nj-.  This 
result  suggests  that  a  reasonable  first  estimate  for  cw  chemical  laser  per¬ 
formance  (other  than  evaluation  of  Nj~)  can  be  obtained  by  assuming  rotational 
equilibrium  and  translational  nonequilibrium,  as  was  done  in  Refs.  6-8.  The 
physical  basis  for  the  latter  approximation  is  as  follows.  CW  chemical  lasers 
operate  at  pressures  of  the  order  of  1  Torr  and  are  inhomogeneous ly  broad¬ 
ened.  The  modification  of  line  shape  (i.e.,  hole  burning),  induced  by  lasing, 
is  more  important  than  lasing  induced  departures  from  rotational  equilibrium. 

The  present  results  may  be  contrasted  with  those  for  a  pulsed  chemical 
laser.  A  pulsed  chemical  laser  generally  operates  at  pressures  of  the  order 
of  one  atmosphere,  in  order  to  achieve  high  energy  density  and  is 
homogeneously  broadened.  As  a  consequence,  its  performance  is  insensitive  to 
translational  nonequilibrium,  and  spectral  output  is  determined  from 
considerations  of  rotational  nonequilibrium. 
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APPENDIX  A.  PARTIAL  LIST  OF  SYMBOLS 


B 

c 

D(  ) 

fj 

gj(v)’  GJ,v 
«v,j(v) 

y1 . 

IS,L’IS,P 


XJ,j 

J 

j 

Ir  I#  u 

cd,Ktr ,Krr 
L(v  -  v') 

NJv,  NJtN,NT 

nv(J,v),nv(J),nv 

Tlfp 

p 

p»  p  »p 

O  V 

p 

Rt  ,Rr 
S 

T 

u 

vy 

X.Xj 

X 

XD 

YJ,j 

y 

cv,j 


parameter  defining  amplifier  saturation,  Eq.  (27) 
speed  of  light  in  vacuum 
Dawson  integral,  Ref.  14 

fraction  of  particles  in  rotational  energy  level,  Eq.  (6) 
normalized  gain,  Eq.  (11) 
gain,  Eq.  (8) 

intensity  for  longitudinal  mode  J,j 
net  intensities 

line  shape  and  power  saturation  intensities,  Eqs.  (12)  and 
(32) 

nondimensional  intensity,  Eq.  (lid) 
rotational  energy  level 
longitudinal  mode  number 
line  shape,  Eqs.  (23)  and  (25) 

deactivation,  translational  and  rotational  relaxation  rates 

Lorentzian  distribution,  Eq.  (8) 

normalized  particle  number  densities,  Eq.  (11) 

particle  number  densities,  Eq.  (1) 

total  number  of  lasing  species,  nj  +  nj 

normalized  output  power,  Eq .  (17) 

velocity  distribution  functions,  Eqs.  (3)-(5) 

pressure 

collisional  rate  ratios,  Eq.  (lie) 

collisional  deactivation  rate  parameter,  Eq.  ( 1 lh ) 

temperature 

streamwise  velocity 

thermal  velocity  in  y  direction 

normalized  frequency,  Eq.  (24) 

streamwise  distance,  Fig.  1 

characteristic  diffusion  distance,  Fig.  1 

homogeneous  width  parameter,  Eq.  (24) 

transverse  distance,  Fig.  1 

energy  per  mole  of  photons 

normalized  streamwise  distance  k^x/u 


^5 


! 


J 


n 

4 


.»*  kl|  . 


a1! Jl* 'Ji*  ^4*  ***•* ' J 


Subscripts 


Superscripts 


normalized  diffusion  distance,  kC{jxD/u 

characteristic  rotational  temperature  parameter,  Eq.  (6) 

wavelength 

frequency 

line  center  frequency 

Doppler  width  [full-width,  half-maximum  (FWHM)] 

homogeneous  width  (FWHM) 

longitudinal  mode  separation,  c/2L 

cross  section  for  stimulated  emisssion,  Eq.  (B-8) 

(1  +  It  t)1/2 


end  of  lasing  region 
start  of  lasing  region 
rotational  level 
longitudinal  mode 
maximum  value 

reference  value  or  rotational  relaxation  value 
vibrational  level 


sum  of  number  densities,  Eq.  (11) 
difference  of  number  densities,  Eq.  (11) 


iSs 

I 

I 

K 


It1' 

S'* 


I 

$ 

M\ 


i 

m 
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APPENDIX  B.  PARAMETER  EVALUATION 


Expressions  needed  to  evaluate  parameters  of  interest  are  noted  below. 


B.l  Dopper  width  [full-width,  half-maximum  (FWHM)] 

1/2 

„  r2kTln2/'* 

Av  “  2  ( - = — J  v  -1 

d  xi  2  o  sec 

Me 
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(831. 6/A )(T/300)1/2(20/M)1/2 


where  X  is  wavelength  in  meters  and  M  is  molecular  weight. 


B.2  Homogeneous  width  (FWHM) 
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where  y  ^  is  2-ncT1'  y^  and  y^  is  tabulated  in  Ref.  7.  Values  of  y^  are  given 


in  Table  III. 


B.3  Gas  kinetic  collision  rate  (HF  +  HF  +  M^) 


where 
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Values  of  a^  are  included  in  Table  III.  For  a  billiard  ball 

model,  y.  ■  a.  and  itAv,  ■  k  ,  . 
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TABLE  III.  Homogeneous  width  and  gas  kinetic  collision 
rate  data  for  HF  +  Mj  +  HF  + 


Mi 

Gas  kinetic  rate 

ai  x  10-11 
(K)  ^^/atm-sec 

Homogeneous  widths 

Rotational  y.  x  10  ^ 

Level  J  1 

(K)  *  ^/atm-sec 

HF 

0.77 

6  +  10 

9.1  +  3.3 

DF 

6+10 

5.7  +  3.0 

h2 

1.95 

5+9 

1.1  +  0.73 

n2 

1.00 

6+8 

1.2  +  0.87 

He 

1.28 

All 

0.16 

Ar 

0.86 

6+8 

0.16  +  0.49 

o2 

0.92 

F2 

0.86 

aRefs.  2  and  6. 
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B.4  Vibrational  deactivation  rate 


We  consider  vibrational  deactivation  of  the  form 

k j 

HF (v  +  1)  +  M  HF(v)  + 


(B-4) 


The  net  deactivation  rate  kcj  is  then 


cd  _  1  din  HF (v  +  1) 

P  “  P  dt 
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Values  of  k  ^ / ( R  T)  are  given  in  Table  IV. 
cd  o 


B.5  Rotational  equilibrium 

For  rotational  equilibrium,  the  fraction  of  particles  with  rotational 
energy  J(J  +  l)kTR  is  given  by 


nv(J)  _  (2J  +  l)exp[-J(J  +  1 >9  R] 

n  '  E ( 2 J  +  l)exp[-J(J  +  1)9  T  ^ 

v  J  R 

-  0  R(  2J  +  1  )exp[-J(J  +  1  )6  R  ]  [  1  +  0(9  R)  ]  (B-6b) 


where  the  summation  is  from  J  =  0  to  J  *  »,  9  *  T  /T,  and  T  =  30.16  K  for 

HF.  The  maximum  value  of  f  and  the  corresponding  value  of  J,  associated  with 

a  fixed  value  of  9  are  denoted  f  and  J  ,  respectively.  Typical  values 

R  J  ,Tn  m 

are  included  in  Table  V.  For  a  fixed  value  of  0  ,  the  quantity  f .  provides 

R  J )  m 

a  convenient  value  for  the  reference  quantity  f  .  When  f  *  f  . 

^  r  r  J,m 


fj  (2J  +  1 )exp [-J (J  +  1 )9  R  ] 

f  J  "  “  (2J  +  1  )exp  [-J  (J  +  1)9_] 

i  to  m  m  k 


and  f j  <  1  for  J  *  Jm« 


B.6  Photon  energy 

Photon  energy,  per  mole,  is 
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TABLE  V.  Maximum  value  of  f.  and  corresponding 

o  1  ~  f  T  £  c  i i  a  ^ 


value  of  J  for  fixed  9  . 

K 


9r 

f  , 

J  ,m 

1/10 

2 

0.2634 

1/20 

3 

0.1889 

1/40 

4 

0.1353 

1/60 

5 

0.1106 

1/80 

6 

0.0957 

1/100 

7 

0.0854 
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where  X  is  in  meters. 


B.7  Gain  coefficient 


expressed 


The  gain  coefficient  for  a  P-branch  transition  v  +  1,J  -  1  -*•  v,  J  can  be 

_ j2 


TAvh  r^-  - 

V  ,  J 


2 

cm 

J-sec 


(B-8) 


where  B(v,J,-l)  is  the  Einstein  coefficient For  HF,  the  latter  can  be 
approximated  by 

B( v , J , -1  )  =  3.79  x  1013  O  +  0.063J)  [l  - 


which  is  believed  to  be  correct  to  within  about  10%  for  1  <  J  <  16  and  v  <  6. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermooechanics ,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy ,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -f ield-of -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radlometry,  mlcrowave/mllilmeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  Infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


